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MISCELLANEOUS. 

73. Proposed by CHARLES E. MTEES, Canton, Ohio (February, 1900. p. 54). 

In an ice cream freezer, cream of a homogeneous character and at the uni- 
form temperature of 60° F. is put into a cylinder having a closed base, and the 
whole put into a freezing mixture so as to subject the base and convex surface 
to a constant temperature of 30° F. Eequired the temperature at any point with- 
in the cream after the expiration of a given time. 

Solution by F. H. SAFFORD, Ph. D., Instructor in Mathematics in the University of Pennsylvania. 
In cylindrical coordinates the equation for the flow of heat is 

D t u=c" [D r *u+— B r ii + \d$ * u + D/ m] 

[Byerly's Fourier Series, p. 1]. The axis of z is to be the axis of the cylinder. 
The planes of the upper and lower bases are respectively 2=0 and z=b, and the 
radius of the cylinder is a. By the nature of the problem, the temperature u at 
any time t is independent of <t>, so that the term above which involves ^ is to be 
omitted. From Fourier's Analytical Theory of Heat, pp. 51, 113, c 2 =K/GD, 
where B is density of solid (cream), € is specific capacity, and K is specific 
conductibility. 

Let #=teniperature of medium surrounding sides and base, ft+»n=intern- 
al temperature when t—0; also let v—m when t=0 (1). 

The differential equation becomes 

D t v=c'- (D/v + — D r v+D z 2 v) (2). 

In order to express the surface conditions due to the surrounding mixture we 
have two more differential equations 

D r ii+h(u — li)=0 when r=a, 
D z ii+h(u—lc)=0 when z—b. 
These become 

D r i)-\-hv=0 when r=a (3), 

D z v-\-hz=0 when z— b (4), 

and are explained in Fourier, p. 94. h is the constant of surface transmission. 

The value of v must satisfy conditions (1) (4). Assume that v—T.R.Z 

where T, B, and Z are functions of t, r, and z, respectively : 

.-.«=» S 2 t CE( " K + pK } .A K J ^ K r).B h sm Pk z (5). 

K = l A=l 

(For method, see Byerly, p. 195). J is a Bessel Function while ,«„ and p K are 
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to be determined by (3) and (4) which give, after v is substituted from (5), 

f-K a<7,(/v a)—ahJ (iJ. K a)—0., (6), 

p K baoBp), b-\-bhsinp K b=0 (7). 

It is shown in Riemann (Partial Differential Equations, §97) that (6) has an in- 
finite number of real positive roots, while tables and a diagram in Gray and 
Mathew's BesseVs Functions make their calculation easy. Fourier, p. 270, states 
the same property for (7). 

We use condition (1) in determing A K and B k , as when t=0 (5) must give 

m=m S 2 A K J (,"k r )B*. sin,o A z. 

It remains to obtain values of A K and B K such that 

2 A K J (X r) =1 , S B k sinp K z = 1 . 
it=i \=i 

From Byerly, pp. 229, 121, we find 
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A * = „,».,... nr,. ^ (8), 



^ = ,f >,^'+^I + l) (1 - C °^ 5) (9) ' 

The same excellent authority gives the means of solving cases in which the ini- 
tial temperature is much more complicated than in the present problem. 
In conclusion, the solution of the given problem is 

B=t+m 2 2 . KU +pk } .A K J {„. K r).B k smp K z } 

the constants being determined by (6), (7), (8), (9). 



